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基于 结构 粒 的 警示 传播 算法 收敛 性 分 析 
牛 进 ， 王 晓 峰 ， 林 青 文 


(北方 民族 大 学 计算 机 科学 与 工程 学 院 , 银川 750021) 


摘 要 : 收敛 性 是 评价 信息 传播 算法 性 能 的 重要 指标 ， 信 息 传 播 算 法 求解 可 满足 性 问题 时 ， 命 题 公 式 的 结构 特征 影 
响 算 法 的 收 么 性 。 具 有 复杂 结构 的 命题 公式 ， 信 息 传 播 算 法 不 总 收 仇 ， 为 了 系统 地 对 此 现象 给 予 理 论 解释 。 借 助 于 
结构 业 的 方法 和 技术 ， 提 出 命题 公式 的 结构 炉 模型 及 其 度量 方法 ,计算 随 机 可 满足 性 实例 的 结构 。 警 示 传 播 算 法 
(WP) 作 为 信息 传播 算法 的 基本 模型 ， 分 析 WP 算法 的 收敛 性 对 于 研究 其 他 信息 传播 算法 的 收 化 性 具有 重要 意义 ,分 
析 了 WP 芷 法 收敛 性 与 结构 蚁 之 间 的 关系 ， 给 出 WP 算法 收 伊 的 判定 条 件 。 通 过 实验 分 析 ， 该 方法 有 效 可 行 。 
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Convergence analysis of warning propagation algorithm based on structural entropy 


Niu Jin, Wang Xiaofeng, Lin Qingwen 
(School of Computer Science 及 Engineering, North Minzu University, Yinchuan 750021, China) 


Abstract: Convergence is an important index for evaluating the performance of information propagation algorithms. When 
the information propagation algorithm solves the satisfiability problem, the structural characteristics of the proposition 
formula affect the convergence of the algorithm, propositional formulas with complex structures do not always converge on 
information propagation algorithms. For this phenomenon, there are few systematic theoretical explanations. As the basic 
model of the information propagation algorithm, the warning propagation algorithm (WP) , analyzing the convergence of the 
WP algorithm is of great significance to study the convergence of other information propagation algorithms. With the help of 
structural entropy methods and techniques, people proposed a structural entropy model of propositional formula and its 
measurement method. This paper calculated the structural entropy of a random satisfiability instance, analyze the relationship 
between the convergence of the WP algorithm and the structural entropy, and give the judgment condition of convergence of 
WP algorithm. Through experimental analysis, this method is effective and feasible. 
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0 ”引言 算法 的 研究 给 出 重要 参考 。WP 算法 是 一 种 基于 因子 图 的 消 
息 迭 代 算 法 ， 通 过 因子 图 的 边 进行 信息 传递 ( 称 为 警示 信息 )， 

约束 可 满足 性 问题 (constraint satisfaction problem) 普 遍 存 ”警示 信息 表示 变 元 的 取 值 是 否 能 够 完全 决定 子 句 的 可 满足 性 ， 

在 于 现实 世界 中 ， 其 研究 目前 在 人 工 智 能 领域 中 具有 广泛 影 。” 该 算法 利用 一 种 迭代 策略 更 新 警示 信息 ， 直 到 算法 收敛 时 停 
响 上 3， 作 为 一 个 重大 子 问 题 ， 布 尔 可 满足 性 判定 问题 ， 止 迭代 ， 或 达到 最 大 迭代 步 数 时 若 算法 仍 不 收敛 ， 则 强迫 迭 
(satisfiability,SAT) 作 为 首先 被 证 明 的 NPC 问题 由 ， 任 何 算法 ” 代 过 程 结束 ,因此 收敛 性 是 WP 算法 最 为 核心 的 性 质 。 目 前 ， 
都 无 法 在 多 项 式 时 间 内 求 得 准确 解 ， 然 而 ， 现 实 世界 中 的 许 ”WP 算法 的 收敛 性 研究 已 经 获得 了 一 些 成 果 , 文献 [10] 证 明了 
多 问题 都 可 以 规约 为 SAT 问题 571。 因此 ，SAT 问题 一 直 在 。”” 当 因子 图 中 只 存在 一 个 环 路 时 , 该 算法 能 够 有 效 收敛 ,得 出 的 
计算 机 科学 界 和 数理 逻辑 界 被 学 者 广泛 关注 。 该 问题 研究 是 。 值 是 变量 边缘 分 布 的 有 效 近似 。 文 献 [11] 基 于 高 斯 模型 给 出 
否 存在 一 组 布尔 变量 ， 赋 值 给 指定 的 命题 公式 ， 能 使 该 公式 。 了 算法 收敛 的 充分 必要 条 件 , 该 条 件 受 限于 两 个 变量 不 能 同 
可 满足 ， 即 使 公式 为 真 。 由 于 任何 一 个 命题 逻辑 公式 都 可 以 ” 属于 一 个 子 句 。 文 献 [12] 基 于 消息 矩阵 分 析 了 算法 的 收敛 性 
转换 为 一 个 等 价 的 合 取 范式 中 ， 因 此 本 文 只 考虑 合 取 范式 。 给 出 了 一 个 充分 必要 条 件 ,但 局 限于 消息 矩阵 为 半 正 定 矩 阵 。 
20 世纪 80 年 代 ， 物 理学 家 提出 了 一 种 基于 因子 图 的 信 文献 [13,14] 分 析 了 基于 高 斯 模型 的 算法 的 收敛 性 ,但 要 求 用 
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息 传 播 算法 ， 该 算法 的 性 能 十 分 依赖 于 因子 图 的 结构 ， 当 因 ” 半 定 规划 来 检查 收敛 性 ,检查 效率 较 低 且 非常 困难 。 文 献 [15] 
子 图 为 单 连通 结构 (任意 两 个 节点 之 间 只 有 一 条 路 径 ) 时 ， 算 分 析 了 求解 均值 的 算法 的 收敛 性 ,要 求 估 计 有 限 尺度 矩阵 的 
法 一 般 能 够 正确 收敛 。 虽然 信息 传播 算法 用 于 求解 SAT 问题 。 谱 半 径 , 在 实际 应 用 中 可 行 性 不 高 。 文 献 [16] 中 分 析 了 算法 的 


时 往往 具有 良好 的 性 能 表现 ， 但 随 着 因子 图 结构 趋 于 复杂 ， 收敛 性 ， 并 给 出 了 信息 传播 算法 收敛 的 一 个 概率 条 件 ， 但 该 
尤其 当 其 中 出 现 环 路 时 ， 信 息 可 能 在 环 路 中 无 限 循环 传播 ， 条 件 主要 局 限于 植 入 指派 的 随机 可 满足 实例 产生 模型 Www” 。 
该 过 程 将 导致 信息 最 终 无 法 收敛 于 一 个 固定 值 ， 因 此 导致 算 。 在 PW” 模型 中 要 求 概率 p 充分 大 ， 其 本 质 是 高 概率 地 最 多 只 
法 无 法 求解 , 该 现象 称 为 算法 不 收敛 四。 WP 算法 作为 一 种 最 个 指派 满足 该 实例 。 不 难看 出 ,上 述 关 于 警示 传播 算法 的 
为 基础 的 信息 传播 算法 ， 研 究 其 收敛 性 可 以 对 其 他 信息 传播 ” 收敛 性 分 析 都 是 基于 某 个 具体 模型 ,或 者 使 用 范围 窄 , 或 者 判 
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定 条 件 的 验证 复杂 等 。 因此 ,对 警示 传播 算法 的 收敛 性 分 析 仍 
然 不 完善 ， 尤 其 基于 普通 因子 图 结构 的 算法 收敛 性 研究 尚 为 


缺失 。 


1.2 WP 算法 


图 的 每 条 边 (4,i) 上 ， 定 义 WP 算 法 中 的 消息 


S 


ai 


警示 信息 ) 为 从 子 句 节点 4a 向 变 元 节点 i 传递 的 一 个 布 


于 WP 算法 的 收敛 ; 


布 的 静态 度量 方法 ,在 信息 论 中 , 若 根据 概率 分 布 选择 节点 ， 
则 焙 值 反映 了 确定 该 节点 的 代码 所 需 藤 平均 信 | 息 量 。 所 
在 通信 和 领域 中 ， 信 息 炉 仅 用 于 度量 点 与 点 之 间 的 单一 信 
在 此 基础 上 , 文献 [18] 首 次 提出 了 结构 丧 。 与 信息 粹 不 同 
结构 炉 是 一 种 动态 度量 方法 ， 可 以 反映 网 络 的 动态 复杂 
因此 ， 结 构 烂 弥补 了 信息 粒 无 法 确切 反映 网 络 结构 复杂 | 
缺陷 ， 可 以 用 以 度量 节点 的 交互 复杂 度 和 内 部 结构 复杂 
重要 信息 。 WP 算法 在 因子 图 上 
杂 性 和 不 确定 性 等 重要 结构 信息 ， 还 可 以 反映 WP 算法 
随机 游 走 过 程 的 动态 复杂 性 。 
综 上 所 述 ， 为 了 研究 WP 算法 求解 SAT 问题 的 收敛 性 ， 


E: 并 全 如 王 


a 
To 
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生 十 分 依赖 于 因子 图 的 结构 ， 而 随 
着 因子 图 规模 的 增 大 ， 因子 图 上 的 节点 之 间 相 互 传递 消息 的 
过 程 不 仅 变 得 极为 复杂 , 同时 具有 随机 不 确定 性 。Shannoni 
提出 的 信息 粒 是 关于 不 确定 性 的 一 种 度量 ， 作 为 关于 概率 分 
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;可 以 表明 变 元 的 取 值 情况 能 否 完全 决定 子 句 4 


wx =1 表 示 : 变 元 5 的 取 值 情况 可 以 决定 子 句 


性 ;uw=0 表示 : 变 元 的 取 值 情况 对 子 句 4 的 
可 满足 性 无 法 起 到 决定 性 作用 ， 即 其 他 变 元 的 取 值 也 能 够 令 
句 a 可 满足 。WP 算 法 的 消息 更 新 迭代 方程 如 下 : 


Ua (t)= II s-4 > JJ it 一 (1) 


jeV(a)i beV (Na 


一 般 改写 为 如 下 形式 : 
a (1)= IT dl 中 Lo (pS 之 ， Up (9)] (2) 


hwa= >》 wt-l)- >》 w(t-1) 03) 


beV*(jN\a beV-(j)\a 


其 中 , hj 为 腔 域 ,车 变 元 尺 只 在 子 句 a 中 出 现 , 则 令 思 =0。 
求解 CNF 公式 F 的 WP 算法 如 下 : 


警示 传播 算法 


本 文 基于 结构 炉 的 概念 ， 提 出 了 命题 公式 的 结构 信息 。 利 用 


G=(m3,m) 模型 随机 生成 不 同 规模 的 SAT 实例 ， 获 取 实 例 的 
结构 信息 。 结 合 WP 算法 在 实例 上 的 收敛 情况 ， 分 析 WP 算 
法 的 收敛 性 和 结构 箭 之 间 的 关系 ， 进 而 利用 结构 烂 度量 WP 


算法 的 收敛 性 ， 给 出 基于 结构 灼 WP 算法 的 收敛 条 件 。 
1 ”基础 知识 


对 于 一 组 给 定 的 布尔 变量 V={w%,%…, 加 } ,集合 中 的 每 一 
个 元 素 称 为 变 元 ， 若 变 元 在 子 句 中 以 自身 值 出 现 ， 对 应 的 
文字 为 *( 正 文字)。 若 % 在 子 句 中 取 反 出 现 ， 对 应 的 文字 为 
ww ( 负 文 字 )。 文 字 的 析 取 可 以 构成 子 句 ， 记 作 
C=vzv…vz)， 每 个 子 句 中 包含 的 所 有 非 重复 文字 个 数 
称 为 子 句 的 长 度 , 记 作 |cl 91。 而 子 句 的 合 取 构成 合 取 范式 公 
式 (CNF)， 记 作 =(G AGA 人 …AG,)， 简 记 为 fF=(G,G,…,G,)。 
SAT 问题 研究 是 否 存在 一 组 对 变 元 的 赋值 指派 ， 使 每 一 个 子 
句 都 为 真 ， 即 使 CNF 公式 为 真 。 特 别 地 ， 当 一 个 CNF 公式 
中 任意 子 句 的 长 度 都 为 , 则 该 CNF 公式 称 为 >- CNF。 由 于 
任何 一 个 CNF 公式 都 可 以 在 多 项 式 时 间 内 规约 为 一 个 3- 


CNF 公式 中 ， 所 以 本 文 只 针对 3-CNF 公式 进行 研究 。 
1.1 因子 图 


设 CNF 公式 为 : 了 =(G,Gw0n) , 变 元 集合 记 为 


fp ， 用 i 表示。 公式 F 可 以 构造 一 个 二 部 图 


ES 


结构 炉 可 以 用 来 度量 WP 算法 在 因子 图 上 的 动态 复杂 性 。 
因子 图 上 的 消息 传播 情 
妈 中 子 句 节点 与 变量 节点 之 间 可 以 通过 边 直接 相连 ， 而 子 名 


G=(CUX,E) ， 称 G 为 因子 图 。 其 中 ,和 = 由 2… 叶 称 为 变 元 集 


合 , C={G,G,…C,} 称 为 子 句 集合 。 
羽 子 图 中 有 两 类 边 : 
实 边 : (C, 力 es 无 全 子 句 C 
虚 边 : (CG, 站 eB 今 子 句 
例如 : 


正文 
全 负 文 


字 上 5 
字 互 


下 =( 人 VDV 一 0) 和 人 (OOV 一 DY 和) 入 


(BVARVXE)AN VAYV XN)= ,043,44 


的 因子 图 G=(V,E) 如 图 1 所 示 。 


Fig.1 Factor diagram 


CNF 公式 F 因子 图 G=(V,E) ， 最 大 迭代 步 tmax 
: if tmax 步 后 WP 算法 未 收敛 ， 
输出 UN-CONVERGED 


输出 所 有 “警示 信息 ”ws 


a) 在 初始 时 刻 t=8, 对 因子 图 G=(V,E) 中 每 一 条 边 , 分 别 以 1/2 的 
岂 从 {96，14} 中 取 一 个 值 赋 给 uw (0) ; 
b) for t=1 to tmax 
顺序 (如 : 字典 序 )， 对 图 G 中 的 每 一 条 边 调 
式 (1) 更 新 ws ; 
(b)if wy(f)=wu s(t-1): 
= Ua (t); 


输出 UN-CONVERGED; 


况 复杂 ， 具 有 随机 性 ， 并 且 由 于 因子 


eo 其 他 边 间 接连 接 ， 同 样 ， 变 量 节 点 2 
间接 相连 。 节 点 之 间 这 样 的 特殊 关系 令 因子 图 产生 
， 因此 ， 本 文 引 入 结构 灶 度 量 此 类 多 层次 因子 
给 定 无 向 连通 图 G=(V,E) ， 设 T 是 G 的 分 割 树 ， 

Laz#4， 结 构 炳 定义 如 下 : 


Hr(G)= 区 (G;w)= i Vv (4) 


是 集合 五 内 部 的 节点 与 该 
a 的 体积 ， 即 7 中 所 有 省 点 的 度数 之 和 。 


V, 


之 外 的 节 ， 


图 G=(V,E) )，8 是 集合 7 
之 外 的 节点 之 间 的 所 有 边 的 权重 之 和 。 图 G 中 节点 vsy 的 度 
由 的 所 有 边 的 权重 之 和 。 


带 权 无 向 连通 图 G 的 K 维 结构 炉 定 义 如 下 : 


H*(G)=min{H™(G)} (5) 

特别 地 ， 带 权 无 向 非 连通 图 G 的 天 维 结构 入 定义 如 下 : 
He (Oit6) Yl(G) Hs (0) GO 
I 是 Ci 的 体积 人 ， 而 Gi， G,, Gi 是 G 的 所 | 
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2 ”命题 公式 的 结构 炳 

根据 WP 算法 收敛 性 分 析 可 知 , 当 因 子 图 不 是 树 型 结构 ， 
而 是 其 中 包含 多 个 环 路 的 一 般 图 状 结构 时 ， 算 法 的 收敛 性 无 
法 保证 ,因此 本 文 利用 高 维 结构 烂 度量 因子 图 的 结构 复杂 性 ， 


从 而 分 析 研 究 WP 算法 在 含有 多 环 路 的 一 般 结 构 因 子 图 上 的 
收敛 性 表现 。 
2.1 无 向 图 转换 技术 


为 了 计算 命题 公式 的 结构 烂 ， 需 要 把 命题 公式 的 因子 图 
模型 转换 为 带 权 值 的 无 向 图 模型 。 因 子 图 利用 两 类 节点 映射 
CNF 公式 中 的 子 句 和 变量 , 并 用 虚 边 和 实 边 表示 节点 之 间 的 
关系 ， 这 样 的 结构 可 以 最 直观 表达 因子 图 和 对 应 公式 的 映射 
以 及 WP 算法 的 运算 过 程 。 构 建 邻接 矩阵 ， 将 因子 图 中 的 两 
类 节点 映射 为 无 向 图 中 的 一 类 节点 ， 以 权 值 代替 因子 图 中 的 
节点 关系 表达 方式 ， 具 体 地 讲 : 

设 变量 节点 为 a 和 b， 子 句 节 点 为 1 和 j， 根 据 式 (1)， 不 
妨 令 ws =1 ， 则 WeV(a)\i， s-[ > wm | 即 


beV(jN\a 


| > jj 于 是 存在 两 种 情况 ， 四 : 7 =-] 且 /=1， 
beV(jN\a 


@: 有 =1 且 =-1。 上 述 两 种 情况 表明 ， 变 量 j 在 子 句 a 和 
中 正 负 出 现 的 情况 相反 ， 根 据 子 句 的 互 异性 ， 在 情形 中 中 ， 
存在 (]) J =1HJ?=1 或 (2) Jr=-1lHJ?=-1 或 者 (3) Jr=1HJ?=-1, 
在 情形 @ 中 ， 存在 (1) J 1 =1 或 Q2) J*=-1H?=-1 或 者 (3) 
Jf=-1 且 7? =1。 上 述 分 析 表 明 ，wy()=1 时 , 环 路 中 正 负 文字 
的 比率 8 存在 两 组 情形 ， 每 组 比率 都 分 别 为 0.25，0.25 以 及 
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7= 人 2… 中 ,节点 按照 地 的 递增 顺序 自 左 向 右 排 序 。 因 此 , 当 


且 仅 当 i<7 ， 将 吃 为 wx 的 左 节点 记 为 ci<Za ; 


gj 对 每 一 个 节点 wsT ,集合 元 cV 。 对 了 
QcB 表 示 4 是 8 的 初始 


的 最 长 初始 字符 段 ; 
d) 对 于 每 个 i， 


字符 段 。 当 节点 gz4， 
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节点 4 和 8B， 


用 w- 表示 


{h(a)= 站 是 V 的 一 种 分 割 情 况 ， 其 中 


A(o) 是 < 的 高 度 ( 根 节点 4 的 高 度 为 0， 并 且 对 于 每 个 节点 


QazA4, h(a)=h(o )+1; 


e) 对 于 每 个 a， 当 


T=Up =QTg; 


P= 时 ， 区 是 7 的 并 集 ， 


人 对 于 T 的 每 个 叶 节 点 & ,是 一 个 单 例 节点 


包含 V 的 单个 节点 。 


因此 ， 


» 大 此 > 1 


初始 化 分 割 树 后 ， 为 计算 结构 焙 ， 对 带 权 无 向 图 进行 最 


小 割 分 割 可 以 得 到 两 个 非 重合 社 
作为 分 割 树 的 第 二 层 树 节点 ,其 


子 图 中 的 所 有 节点 ， 


M 包含 所 节点 的 数 
进行 分 割 ， 


M, = M, J M,, 。 若 子 社 


即 所 有 节点 划分 为 一 个 社 
一 个 社区 可 能 包含 一 个 或 多 个 节点 ， 且 V=M UM,。 
大 于 1， 则 继续 利 


继 节 点 ，M 和 M> 为 M, 的 直接 后 


区 M 和 M,， 将 这 两 个 社区 
父 节点 都 为 根 节点 


4 (包含 因 


又 )， 


中 ,每 
若 妇 或 


最 小 割 对 2 或 4 
将 分 割 可 能 得 到 的 子 社区 Mu ，M 以 及 Ma ，M> 
作为 分 割 树 的 第 三 层 树 节点 , 其 中 ，WM 和 Mi; 为 Wi 的 直接 后 


全 节点 ， 且 Mi=Mi JIM, ， 


区 中 仍 包含 一 个 或 多 个 节点 , 则 重复 如 


上 步骤 继续 分 割 子 社区 ， 并 构造 分 割 树 ， 直 到 当 社 区 之 中 只 
包含 单个 节点 时 ， 对 该 社区 停止 分 割 。 
若 因子 图 为 带 权 无 向 非 连 通 图 ， 则 遍历 其 节点 ， 找 到 每 


个 连通 分 量 ， 对 其 每 个 连通 分 量 进行 如 上 的 图 分 割 并 且 构 


造 分 割 树 。 


0.5。 大 wx*(O)=0， 同 理 可 证 ， | > Wm ]>0 时 ， 存在 两 
beV(jNa 


种 情况 ，@: 下 =-1 且 六 =- ， 加 ， 好 =1 且 及 =1。 根 据 上 述 分 
析 ，B 存在 两 组 情形 ， 每 组 比率 都 分 别 为 0.25 以 及 0.25。 而 


当 必 | 郧 /]-0 时 ， 由 于 结果 只 受 %1 影 响 ,因此 包 合 
beV (iNa 


以 上 所 有 可 能 性 ， 且 对 最 终结 果 不 会 产生 影响 ， 所 以 不 对 这 
种 情况 做 具体 分 析 
显然 ， 在 基本 环 路 模型 中 ， 正 文字 与 负 文 字 的 比率 期 望 


值 为 E(oO=03 ， 将 该 期 望 值 设 为 负 文字 或 正文 字 的 权 值 。1 
于 结构 粒 只 关心 因子 图 的 结构 信息 , 无 论 虚 边 权 值 为 0.3、 实 
边 权 值 为 0.7, 还 是 采取 相反 赋值 ,两 种 赋值 产生 的 图 仅 为 对 


称 结构 ， 不 影响 结构 信息 。 因 此 ， 本 文 为 了 统一 结果 ， 均 设 
虚 边 权 值 为 0.3， 实 边 权 值 为 0.7。 采 用 此 方法 ， 将 图 1 因子 
图 转换 为 带 权 无 向 图 如 图 2 所 示 。 


图 2 带 权 无 向 图 
Fig.2 Weighted undirected graph 

结构 类 模 型 及 其 度量 方法 
给 定 由 7 个 节点 和 mw 个 边 组 成 的 加 权 无 向 连通 图 
G=(V,E) ， 为 了 计算 结构 米 ， 利 用 社区 划分 构造 由 图 中 节点 
集合 作为 树 节点 的 分 割 树 ， 该 树 型 结构 可 以 根据 以 下 规定 进 
行 构造 ， 其 中 ，P?={Xi,X…Xi} 是 节点 集合 Y 的 分 割 ，X; 称 
为 社区 : 

a) 将 根 节点 表示 为 +4， 定义 节点 集合 T=V ; 

b) 对 于 每 一 个 节点 geT ,2 的 直接 后 继 节 点 为 a"j， 


2.2 


因子 图 的 最 小 割 过 程 与 对 应 分 割 树 


图 3 和 4 所 示 ( 为 方便 表示 ， 仍 然 


的 构造 过 程 如 
虚线 和 实 线 蔡 代 权 值 ):; 


on © 
XO 
GY GY 


中 


(9 


© © OO OO 


@ 


区] 


3 


刀子 图 分 割 


Fig. 3 Factor graph segmentation 


(7 ) T={1,2,3,4,5} 


4 


分 割 树 构造 


Fig.4 Partition tree Structure 


分 割 树 构造 完成 后 ， 利 
中 节点 的 度数 为 对 应 带 权 无 向 图 
节点 的 米 值 作 和 得 到 #7(G) 。 
炉 应 取 87(G) 的 最 小 值 ， 即 找到 一 种 最 优 分 割 。 


的 焙 值 ， 其 
数 ， 对 所 


] 式 (4) 计 入 


分 割 树 中 的 每 个 节点 
FPF 相应 节点 的 度 


前 学 术 界 


202009.00095v1 
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对 于 社区 发 现 的 研究 正在 进行 ， 由 了 
区 发 现 ， 社 区 之 间 


大 ， 所 以 最 小 割 自 
泛 己 ]。 利 用 在 前 节 


的 边 密度 较 小 ， 社 区 内 部 的 边 密度 本 
法 在 现 阶段 入 


ml 


[ll 


要 影响 ,所 以 本 实验 利 


值 对 边 赋予 权 值 ， 


中 求 得 正 负 文字 比率 的 期 
比 ， 基 于 该 权 值 的 最 小 


满足 结构 焙 
中 的 结构 信息 损耗 。 
即 为 CNF 公式 的 结构 信息 

算法 2 CNF 公式 结构 信息 获取 


输入 : 
输出 : 


begin 
a) XO{2,,n)}, COO{O,G,.…C,} ， 构 j 


的 社区 定义 ， 


CNF 公式 F 
H*(G) 


G=(CUX,E); 


b) 实 线 FL=8.7, 虚线 DL=8.3， 构 造 带 权 无 向 
c) 初始 化 分 割 树 ， 了 到 = 纪 ; 

d) if G(F) connect: 

go to step e; 


flag = 1 


else: 


e) 利 


遍历 所 有 连通 分 量 ，go to step ej 


flag = 6 


go to step 8g; 


else: 


go to step el; 


8g) 利 


最 小 割 结果 构造 分 割 树 ; 


h) if flag = 1: 


H* (G)=min{H™(G)} , 


else: 


end 


再 加 


播 过 程 ， 并 且 ， 基 于 该 过 程 所 产生 的 社 


H*(G,)= > HT™(G;a) 


aeT,az#A 


天 十 1 总 D .HK 
H (0 = yo) lls H*(G,) 


i) return H*(G); 


上 述 方法 所 求 结果 为 wr (6)}， 


变 元 个 数 分 别 设置 为 n=10 和 
数 为 100。 以 n=10 的 情况 为 例 ， 
定 4 的 规模 ， 利 用 模型 控 人 
随机 实例 。 为 使 实验 结果 更 加 精准 
生 的 实例 个 数 都 为 100， 并 对 每 一 个 实例 的 因子 图 进行 上 文 
所 述 分 割 ， 最 终 根据 算法 2 计算 得 到 该 命题 公式 的 结构 和 。 


通过 固 


: 基于 结构 蚁 的 警示 传播 算法 收 伊 性 分 析 


WP 算法 的 收敛 性 分 析 
于 比值 a=m/n 对 因子 图 


厂 
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的 结构 和 WP 算法 收敛 性 具 


其 中 ,代表 变 元 集合 中 所 有 变 元 
k=3 代表 每 个 子 句 都 恰好 包含 3 个 变 元 。 实 验 选 取 两 组 数据 ， 


模型 G=(m,3,m) 产生 随机 实例 。 
的 个 数 ,到 代表 子 句 个 数 ， 


n=20 


， 令 WP 算法 最 大 迭代 次 


体 实验 过 程 为 : 令 k=3， 


%， 用 以 产生 不 同 规模 的 


， 其 中 ， 令 不 同 规模 下 产 


了 计算 100 个 实例 的 结构 粒 的 均值 作为 该 实例 组 的 结构 箭 。 
随 着 a 从 0 增 至 5， 随 机 产 4 


最 小 割 分 割 图 ; 
f) if Subgraph == single node: 


日 式 (4) 可 知 ， 基 了 
业 值 相 加 所 得 ， 即 ; 
子 图 的 结构 炉 考 虑 了 WP 算法 在 


以 有 诡 


度量 社 


区 之 间 的 动态 信息 交互 和 社区 内 部 的 信息 交互 


过 程 ， 


而 子 社区 的 分 割 过 


区 分 嘎 


定性 等 


| 情况 ， 


象 为 结构 灶 的 过 程 中 的 信息 损耗 , 该 算法 求 得 CNF 公式 结构 


厅 了 因子 图 的 结构 复杂 性 。 结 构 和 
的 动态 复杂 性 和 不 确 


从 根本 | 
可 以 有 效 并 完整 地 分 析 者 
EE 要 性 质 。 利 用 
最 小 制 ， 可 以 有 效 去 除 图 


司 时 最 小 化 结构 信息 


炉 , 其 中 包含 着 CNF 公式 的 结构 复杂 性 , 不 确定 怕 
的 计算 结果 是 CNF 公式 的 结构 信息 
并 以 结构 灶 


将 该 算法 用 于 分 析 WP 算法 收敛 性 研究 中 ， 不 仅 不 用 霸 
模型 ,使 用 
以 往 WP 算法 收敛 和 
因子 图 结构 的 WP 算法 收敛 性 丰 


性 等 重要 结构 信息 , 最 终 
的 高 度 抽象 ，j 

于 特定 结构 的 图 

并 且 可 以 填 率 


WP 算法 
炉 以 实验 验 订 


研究 提供 


新 的 思路 与 方向 。 接 下 来 本 文 将 霸 
F 分 析 WP 算法 收敛 性 。 


E 的 3SAT 实例 的 结构 业 变 


化 情况 如 图 5 和 6 所 示 。 其 中 ， 每 个 数据 点 的 结构 炉 均 由 每 


4 


网 组 中 所 有 实例 的 结构 炉 
含 100 个 随机 实例 。 
算法 成 功 收敛 的 实例 
的 概率 随 结构 灼 变化 如 菇 


同时 ， 太 
个 数 ， 得 


计算 均值 得 出 ， 每 个 实例 组 包 
[于 每 个 实例 组 ， 统 计 其 中 WP 
1 该 实例 组 的 收敛 概率 。 实 例 
7 和 8 所 示 。 
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图 5 结构 炉 分 布 (n=10) 
Fig.5 Structural entropy distribution( n=10) 
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图 6 结构 炳 分 布 (n=20) 


Fig.6 Structural entropy distribution( n=20) 
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图 7 收敛 性 分 布 (n=10) 


Fig.7 Convergent distribution( n=10) 


ChinaXiv 合 作 期 刊 
录用 定稿 牛 进 ， 等 : 基于 结构 炳 的 警示 传播 算法 收 伊 性 分 析 第 38 卷 第 3 期 
嘱 A Pe constraint satisfaction problems [J]. Discrete Mathematics, 2009, 309 (8): 
四 2085-2099. 
| i [3] Martin OC,Monasson R,Zecchina R. Statistical mechanics methods and 
3 "| phase transitions in optimization [J]. Theory Comput Sci, 2001, 265: 3— 
-| ! 67. 
? | T [4] Cook, Stephen. The Complexity of Theorem-Proving Procedures [C1// In 
* Proceedings of the third annual ACM symposium on Theory of 
2 | T computing, 1971: 151-158. 
oz| i [5] Clarke E, Biere A, Raimi R, et al. Bounded model checking using 
6 ; 二 4 二 ， 志 satisfiability solving [J]. Formal Methods in System Design, 2001, 19 (1): 
Structure Entropy(H) 7-34. 
图 8” 收敛 性 分 布 (n=20) [6] Vizel Y, Weissenbacher G, Malik S. Boolean Satisfiability Solvers and 
Fig.8 Convergent distribution (n=20) Their Applications in Model Checking [J]. Proceedings of the IEEE, 
实验 表明 , 随 着 4 增 大 , 随机 实例 的 结构 粒 也 随 之 增 大 ， 2015, 103 (11): 2021-2035. 
但 结构 灶 增 量 随 之 减 小 。 图 $ 可 知 ， 当 0<w<2 时 ， 结 构 粒 [7] Kuper L, Katz G, Gottschlich J, et al. Toward Scalable Verification for 
增长 较 快 ， 当 w<>2 时 结构 业 增 幅 较 为 平缓 ， 特 别 地 ， 当 cx 趋 Safety-Critical Deep Networks [J]. 2018, 34 (5): 312-322. 
近 5 时 ， 结 构 和 增 量 趋 近 于 0。 图 7 表明 <w>3 时 ， 结 构 灶 变 [8] Braunstein A, Mezard M, Zecchina R. Survey propagation: an algorithm 
化 较为 平 组 6 for satisfiability [J]. Random Structures & Algorithms, 2002, 27 (2): 201- 
当 n=10 时 ， 随 着 “ 趋 近 于 5， 实例 的 炉 值 密集 分 布 于 226. 
3.5324-3.8941 之 间 , 分 布 在 该 区 间 的 实例 占 实例 总 数 的 74%， [9] Weiss Y. Correctness of Local Probability Propagation in Graphical 
最 终 炉 值 允 近 3.9。 当 业 值 小 于 3.7148 时 ,WP 算法 在 所 有 实 Models with Loops [J]. Neural computation, 2000, 12 (1): 1-41. 
例 上 完全 收 人 鳅 ， 炉 值 大 于 3.7148 时 ， 算法 收敛 概率 从 完全 收 [10] Weiss Y, Freeman W T. Correctness of Belief Propagation in Gaussian 
敛 减 小 到 0.71。 Graphical Models of Arbitrary Topology [J]. Neural Computation, 2001, 
当 n=20 时 , 实例 的 灶 值 密集 分 布 于 3.8136-4.3494 之 间 ， 13 (10): 2173-2200. 

最 终 逼 近 4.35。 当 焙 值 小 于 4.2795 时 ，WP 算法 几乎 完全 收 ” [11] 王晓峰 ， 许 道 云 ， 书 立 , 随机 可 满足 实例 集 上 警示 传播 算法 的 收敛 
敛 ， 当 焙 值 大 于 4.2795 时 ， 算 法 收敛 概率 急剧 下 降 ， 最 终 从 性 [J]. 软件 学 报 , 2013 (01): 5-15. (Wang Xiaofeng, Xu Daoyun, Wei 
1 减 小 至 0.66。 Li. Randomly satisfy the convergence of the warning propagation 

结构 灶 的 变化 情况 表明 » 当 变 元 规模 国定 时 ， 随 厦 wx 的 algorithm on the instance set [J]. Journal of software, 2013 (01): 5-15) 
首 大 ， 因子 图 的 动态 复杂 性 会 在 一 定 范围 内 快速 上 升 ， 超 过 [12] 王晓峰 ， 许 道 云 ,警示 传播 算法 收 化 的 充分 条 件 [J]. 软件 学 报 ， 

值 时 ， 动 态 复杂 性 的 上 升 速度 会 显著 减缓 。 这 是 由 于 在 简 2016 (12): 3003-3013. (Wang Xiaofeng, Xu Daoyun. Sufficient 
单 因 子 图 结构 中 ， 任 何 的 变化 ， 如 环 路 增多 ， 图 规模 的 细微 conditions for convergence of warning propagation algorithm [J]. Journal 

化 ， 以 及 因子 图 直径 的 改变 ， 都 会 对 因子 图 结构 造成 较 大 of software, 2016 (12): 3003-3013.) 

响 。 而 当 因 子 图 规模 增 大 到 一 定 程 度 时 ， 环 路 的 增加 以 及 [13] Su Qinliang, Wu Y C. Convergence analysis of the variance in gaussian 

子 图 在 一 定 范 围 内 的 规模 增 大 ， 对 于 图 的 动态 复杂 性 影响 belief propagation [J]. IEEE Transactions on Signal Processing, 2014, 62 

j 显 减 小 ， 即 因子 图 状态 更 为 稳定 。 由 算法 的 收敛 性 变化 可 (19): 5119-5131, 

， 当 因子 图 状态 越 为 稳定 时 ，WP 算法 的 收敛 性 随 着 因子 [14] Su Qinliang, Wu Y C. On Convergence Conditions of Gaussian Belief 

到 的 动态 复杂 度 升 高 而 降低 。 而 当 因子 图 的 动态 复杂 性 未 突 Propagation [J]. IEEE Transactions on Signal Processing, 2015, 63 (5): 


破 闵 值 时 ，WP 算法 可 以 高 概率 收敛 。 

妹子 图 的 结构 决定 了 WP 的 性 能 ， 实验 给 出 了 不 同 
规模 下 ，WP 算法 收敛 的 充分 条 件 ， we 若 
H*(G)<3.7148 ，WP 算法 高 概率 收 你 。 当 ”=20 ， 和 若 
H*(G)<4.2795 ,WP 算法 高 概率 收敛 。 基 于 其 他 规模 因子 图 的 


WP 算法 收敛 性 判定 条 件 ， 同 样 可 以 利用 该 方法 得 到 。 
4 ”结束 语 

本 文 提出 了 CNF 公式 结构 信息 获取 算法 , 并 利用 结构 炉 
表达 结构 信息 ， 通 过 实验 得 到 了 不 同 因子 图 规模 下 结构 灼 
变化 情况 以 及 在 不 同 结构 粹 范围 内 WP 算法 在 因子 图 上 的 收 
敛 性 ， 给 出 了 不 同 规模 下 WP 算法 收敛 的 充分 条 件 。 下 一 步 
工作 将 寻找 关于 命题 公式 结构 粒 或 控制 参数 ， 用 以 
度量 其 他 信息 传播 算法 的 收敛 性 。 并 对 结构 焙 和 因子 图 结构 
言 息 之 间 的 关系 做 更 进一步 的 理论 分 析 和 实验 验证 。 
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